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Practically all of our knowledge of the Golgi apparatus has been built 
up through the study of fixed, osmic or silver, preparations, it not being 
generally realized that this important cytoplasmic inclusion can be re- 
vealed in many living cells by means of .vital methylene blue staining. 
Reliance on the fixed preparations has created a great deal of confusion 
with regard to the structure of the Golgi bodies, because their nature is 
such that they only rarely lend themselves to faithful preservation by 
fixing fluids. In molluscs, these inclusions are of such large size that it has 
been possible to follow their behavior in intact embryos over periods of 
several days and considerable new information has been gathered regarding 
their structure and function. 

Such observations reveal that the generalized Golgi body in molluscs 
is typically a more or less spherical, simple vesicle, possessing a chromo- 
philic, gel-like, relatively lipoidal pellicle covering a relatively more pro- 
tein, fluid, chromophobic core. The pellicle is thickened over one surface 
of the sphere forming a bowl-shaped structure which, through the micro- 
scope, looks like a crescent (Fig. 1 (A)). The relation between the thick- 
ened, bowl-shaped portion of the chromophile and the chromophobe is 
somewhat like the relation that exists between the “crescent” moon 
and the moon as a whole, except that in the Golgi body a thin layer of 
chromophilic material completely encloses the chromophobe. These gen- 
eralized bodies are clearly endowed with the ability to keep the lipoidal 
and protein components segregated, hence the visible differentiation into 
chromophile and chromophobe. But they are the immediate descendants 
of homogeneous droplet-like forms where no such ability is manifest (Fig. 
1 (B)). These droplets, in turn, are derived from smaller and presum- 
ably more solid granules which may be considered ‘‘reserve” forms (Fig. 


1 (C)). 
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The enlargement of the granules to form droplets and the development 
of the droplets into larger, simple vesicles are due, with very little doubt, 
to the absorption of substances from the cytoplasm. This absorbing 
quality, together with the segregating tendency already referred to, con- 
stitute the two most fundamental physiological characteristics of the Golgi 
substance in these species. 

In developing molluscs, the generalized Golgi bodies give rise, following 
gastrulation, to exceedingly active and highly specialized forms. By this 
period, quantities of the original yolk have been dissolved in the cyto- 
plasm and although much has, no doubt, been oxidized, the Golgi elements 
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FIGURE 1 
Living Golgi bodies seen in trochophore of the tectibranch Navanax 
(X 1000). A, generalized form, showing “crescent’’-like chromophile, 


and chromophobe; B, droplet; C, granule; D, compound vesicular 
forms, each showing several chromophobes. 


absorb large amounts of the material and increase considerably in size, 
forming elaborate, compound vesicular bodies (Fig. 1 (D)). Within the 
chromophilic material, new vesicles continually make their appearance 
and, in the trochophore larva, small (1”) fat droplets are continually dis- 
charged into the cytoplasm. In the pelecypod Mytilus, the absorption by 
the Golgi bodies of all of the ultramicroscopic, pigmented fat particles of 
the single-celled egg results in the freeing of the ground cytoplasm of pig- 
ment and the deposition of this pigment in the Golgi inclusions. 

The discharge from each enlarged Golgi body of large (4-6”) protein 
spheres (the original chromophobes) is characteristic, in nudibranchs and 
tectibranchs, of the early veliger larva. In this process, the chromophilic 
portion of the structure withdraws from around the chromophobe and 
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condenses into a chromophilic, pycnotic mass, which can be shown to be 
rather solid. The pycnotic phase, however, is soon followed by a stage in 
which the mass becomes resolved into a number of spherical, relatively 
fluid, droplets. Each of these soon differentiates into a generalized vesicle 
as the segregating tendency is resumed and the cycle of the system begins 
anew. 

The breaking up of the pycnotic Golgi body into droplets is the repro- 
ductive phase of the cycle and is the only method by means of which the 
Golgi elements increase in number during early development. In other 
words, there is no “‘dictyokinesis” accompanying every division of the cells, 
all multiplication being a mass division of the sort just described and al- 
Ways occurring after the assumption of a pycnotic condition following a 
period of marked synthetic activity. This also appears to be the situation 
during odgenesis where the synthetic period results in the deposition of 
large quantities of fat and protein yolk in the oécyte. The multiplication 
period clearly increases the number of Golgi elements, thereby in tecti- 
branchs providing a sufficient number for the new cells of the growing 
embryo. In the event that the embryo is not increasing in size during this 
period, as in Mytilus, this increase in number may be offset by a fusion 
phase preceding the active synthetic period. The fusion phase reduces the 
total number of Golgi bodies, but the original number is approximately 
restored during the multiplication period. 

In the account just given, the internal chromophobe of each Golgi 
body is a potential protein or fat product which gradually becomes much 
more viscous as development proceeds. The elaboration of such products, 
in molluscs, occurs chiefly within the Golgi vesicle, but elaboration outside 
the vesicle is employed in Mytilus in the development of certain large 
oily yolk droplets. In this case, the product first makes its appearance 
within a group of small Golgi bodies, all of which appear to contribute to 
the single common yolk sphere. The mechanism involved in this process 
is not clear, but it is significant that elaboration of a product can take 
place merely through contact of the Golgi body with the ground cyto- 
plasm. 

In developing molluscs, the function of the Golgi system appears to be 
the continual mobilization of the fat and protein reserves within the cell. 
Why the embryo should continually dissolve its formed yolk inclusions, 
only to elaborate new ones is not immediately evident, but that the phe- 
nomenon is of fairly general occurrence in animals is indicated by the work 
of Schoenheimer! who finds that almost all of the proteins of the body are 
continually undergoing synthesis and breakdown. Should the Golgi 
system be found to be concerned in this protein turnover in higher animals, 
as it seems to be in the molluscs, it would prove to be a cytoplasmic con- 
stituent of greater metabolic significance than we have realized. 








228 ZOOLOGY: L. G. WORLEY Proc. N. A. S. 


Summary.—For the first time, the Golgi system has been continuously 
followed in living, developing animal eggs. The multiplication of the 
Golgi elements has been observed and the réle of the Golgi substance as an 
absorbing mechanism and protein and fat elaborating system is described. 
The active Golgi bodies are found to originate from minute, chromo- 
philic ‘‘reserve”’ granules. 


* Contribution from the Department of Biology, Brooklyn College, No. 45, and 
Contribution from The Scripps Institution of Oceanography, New Series No. 208. 
1 Schoenheimer, R., The Dynamic Body Constituents, Cambridge, Mass., 1942. 
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There has been considerable difference of opinion concerning the na- 
ture of certain ‘droplets’ that become visible in cells when the tissue is 
stained with various basic vital dyes. One group of workers, notably 
Beams,! Ludford,*? Chang* and others, holds that these are ‘‘neoforma- 
tions’ resulting from the vital dye treatment and that they have no counter- 
part in the unaltered, living cell. Another group of investigators, notably 
Covell and Scott,4 Cowdry and Scott,’ Ma,* Owens and Bensley’ and Parat,® 
are of the opinion that these “dye droplets’ represent the living Golgi 
material and that the classical Golgi apparatus of the fixed and osmicated 
cell results either from a precipitation of the osmic salt in and around these 
droplets, or from their running together and fusion as the result of fixation. 
In recent years, the tendency has swung in favor of the first contention, 
but crucial proof of neither of these views is at hand. 

In my work, “droplets” of this kind have been revealed by means of 
supra- or intravital methylene blue staining of the tissues of many inverte- 
brates, including the salivary gland cells of the larva of the midge Chirono- 
mus, and the smooth muscle, pancreas, liver and gall bladder of the frog 
and the smooth muscle, adipose tissue, pancreas, liver, kidney, thymus, 
bladder and pars distalis of the kitten. Usually, these ‘‘droplets” are 
opaque and show, in life, no internal structure and the reasons for their 
being considered “‘dye droplets’’ by so many workers are evident. 

However, in all cases thus far examined, when living tissues containing 
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these ‘‘dye droplets” are treated with hypertonic salt solution, thus simu- 
lating the shrinkage that so frequently attends osmic acid fixation, the 
“droplets” rupture and form in each cell a group of fixtures that collec- 
tively so closely resemble the classical Golgi apparatus that there is little 
doubt but that the two are identical. 

Experiments with hypertonic salt solution have shown that each ‘‘excre- 
tion droplet” is in reality a duplex structure, consisting of a rather thick, 








pe Sees 28 


FIGURE 1 
Living Golgi vesicles in the salivary gland of the larva of Chironomus, supravitally 
stained with methylene blue (X 500). 
FIGURE 2 
Same five minutes later after treatment with hypertonic salt solution (X 1000). The 
vesicles have ruptured, producing the classical Golgi figures. Compare with figures of 
Beams and Goldsmith.® 


chromophilic, gel-like pellicle and a fluid core, the latter corresponding to 
the chromophobe or osmiophobe of many invertebrate Golgi bodies. 
The fluid interior appears to be held under some slight pressure by the 
pellicle and the structure, far from being a mere droplet of dye is actually 
a simple turgid vesicle. 

When these Golgi vesicles in the relatively inactive gland cells of Chi- 
ronomus (Fig. 1), for example, are treated with hypertonic salt solution. 
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many chromophilic rings are immediately formed, each of which usually 
shows a localized thickening (Fig. 2). These vesicles, therefore, have an 
internal structure which no ordinary dye droplet could possibly be claimed 
to possess. Furthermore, the pellicle is revealed by these experiments 
to be much thicker than would be the case of a surface membrane that 
might conceivably be accredited to an excretion droplet. When relatively 
active gland cells are similarly treated, the fluid vesicle content squirts 
out into the cytoplasm forming chromophilic ‘‘streamers’”’ which immedi- 
ately solidify, or the entire vesicle may elongate forming a coarse thread. 
Under favorable circumstances, these threads and streamers, in the mam- 
malian acinar cells of the pancreas, contact each other forming a chromo- 
philic pseudoreticulum which bears a striking resemblance to the classical 
Golgi apparatus of the osmicated cell. This resemblance may be further 
increased by allowing the material to desiccate on the slide, thus simulating 
tissue dehydration. Frequently, both intact, and ruptured or elongated, 
vesicles may occur side by side in the same cell. This probably accounts 
for the widespread belief that since “dye droplets’ or ‘‘vacuome (?)” 
(unruptured Golgi vesicles) and Golgi apparatus (ruptured or elongated 
Golgi vesicles) can be demonstrated in the same cell, the two are not iden- 
tical. Indeed, this may be true in some tissues, but it is certainly not the 
case in those mentioned above. 

When vitally stained amphibian smooth muscle cells are treated with 
hypertonic salt solution, the Golgi vesicles become flattened into discs 
which, upon desiccation, due to their soft gel-like consistency often fuse to 
form a longitudinal plate-like structure. It is probably in this manner 
that the classical lamellar form of the fixed Golgi apparatus in these cells 
is brought about by the fixing fluid. When vitally stained tissue is treated 
with hypertonic salt solution after a considerable delay following its re- 
moval from the animal, distortion of the vesicles is less marked, probably 
because on standing the vesicles undergo a considerable increase in viscos- 
ity. 

These statements are not intended to imply that true excretion droplets 
never occur in vitally stained tissues, but it is perfectly clear that in the 
tissues thus far examined, the classical Golgi apparatus of the fixed cell is 
formed from chromophilic vesicles of the living cell that have frequently 
been considered “‘neoformations.”’ It appears very probable that there is a 
correspondence between the ‘“‘vacuome”’ and the Golgi system in far more 
instances than’ we have been accustomed to believe. At the same time, 
the vital methylene blue staining technique has revealed true Golgi pud- 
dingstone-like masses in invertebrate. odcytes and developing embryos 
which, in section, appear network-like. 

Summary.—The living Golgi system in many invertebrate and verte- 
brate cells consists of a series of chromophilic vesicles, each with a rather 
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thick gel-like pellicle enclosing a fluid core. Upon treatment of the vitally 
stained tissue with hypertonic salt solution, thus simulating fixation 
shrinkage, these vesicles can be observed to collapse, discharge the fluid 
vesicle content or elongate forming a pseudoreticulum similar to that ex- 
hibited by the Golgi apparatus of the osmicated cell. Desiccation, super- 
imposed upon the shrunken material, provides an even more striking pic- 
ture of the classical Golgi apparatus. 


* Contribution from the Department of Biology, Brooklyn College, No. 46, and 
Contribution from The Scripps Institution of Oceanography, New Series No. 209. 
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MATTER, ELECTRICITY AND GRAVITATION IN FLAT 
SPACE-TIME 


By Greorce D. BIRKHOFF 
HARVARD UNIVERSITY . 
Communicated July 13, 1943 


In 1927 I presented two Notes! in which there was attempted a con- 
ceptual approach to the then new Schrédinger wave equation. This was 
done by taking matter to be a “perfect fluid,’ defined against the back- 
ground of the curved space-time of Einstein’s celebrated gravitational 
theory of 1916. 

In February of last year I had the honor of presenting at Puebla and 
Tonantzintla, Mexico, before the Astrophysical Congress convening there, 
a new gravitational theory based on the same perfect fluid in the much 
simpler flat space-time characteristic of modern electromagnetic theory 
and special relativity. Already in 1912 Nordstrém had proposed a gravi- 
tational theory founded upon this type of space-time;? but his theory 
failed in that it did not account for the observed slight advance of the peri- 
helion of Mercury beyond the amount indicated by the classical Newtonian 
theory. 

My theory was obtained by laying down those demands which seemed 
most natural and elegant from the mathematical point of view. 
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In the present Note.I wish to indicate in outline the modified account of 
matter, electricity and gravitation thus arrived at. The appropriate 
mathematical language is no longer that of tensors as in my two Notes of 
1927, but is that of 4-vectors. It should be emphasized that, from the 
mathematical and philosophical point of view, the new theory is very 
simple. 


I. NORMAL COORDINATES IN FLAT SPACE-TIME 
Let ds denote the element of local time so that 
ds? = dt? — dx? — dy — dz? 


where dt and dx, dy, dz refer to the usual time and space codrdinates in 
seconds and light-seconds, respectively. If now we replace ¢, x, y, 2 by 


si=t,<?= V —1x, ag V—1y, xt=V—-— lg, this formula takes the form 
ds? = (dx)? + (dx?)? + (dx*)? + dx*)?, 


and the corresponding coérdinates x’ are called normal codrdinates. In 
such normal coérdinates the language of 4-vectors becomes the same as 
that of 3-vectors in ordinary space. For this reason it is possible to use 
subscripts throughout rather than the subscripts and superscripts charac- 
teristic of tensor theory. Thus we have 


ds* = dxa’ (1) 
For brevity we shall use such normal coérdinates almost exclusively. 


II. THe PerFect FLuIp 


By general consent the energy tensor of the homogeneous adiabatic 
fluid is written in essentially the form 





Ty = pugs — pbs (2) 
and, in terms of it, the equations of motion are written 
OT ia 
= fy. 3 
wihetL (3) 


Here p and p = f(p) designate, respectively, the density and pressure of the 
fluid, u; is the velocity vector dx,/ds, 5,; is the usual Kronecker 6,;, and f; 
stands for the body force vector per unit of volume. 

The perfect fluid is singled out by the further requirement that the dis- 
turbance velocity is to be that of light at all densities, with the correspond- 
ing equation of state, 


b = p/2. (4) 
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Only with this type of fluid can essential mathematical difficulties be 
avoided at collision of portions of the fluid.* It is also assumed that there 
is an equilibrium density po,corresponding to a cosmic pressure po/2. The 
precise value of this density nowhere enters, and it would be equally pos- 
sible to suppose that the equation of state has the form » = (p — po)/2 
and that the pressure is 0 at the free boundaries in the customary manner. 

We recall that the force vector f; is necessarily orthogonal to the velocity 
vector u, 


fata = 0. (5) 


III. THe GENERAL HYPOTHESES ON THE FORCES 


We propose now to make the following assumptions: (a) the force vector 
fi is rational and integral in the velocity components, of not higher than the 
second degree; and (bd) the coefficients are homogeneous and linear in the 
first partial derivatives of the corresponding potentials, namely, the atomic 
potential y which I introduced in 1927 (loc. cit.), the usual vector electro- 
magnetic potential y;, and the symmetric gravitational tensor potential 
hy defined in the present Note. We shall furthermore suppose that (c) 
there are no degenerate quadratic terms in the velocities, i.e., no terms re- 
ducing to terms independent of u, in view of the fundamental identity 
ua? = 1, These three types of potential seem appropriately designated, in- 
asmuch as they refer primarily to matter, to electricity and to gravitation, 
respectively. 


IV. THE ATOMIC POTENTIAL y 


The components of the external forces acting on unit volume of the 
perfect fluid due to the atomic potential will, according to our general hy- 
pothesis above, be given by a vector of the form 


“~* OXa 


of degree 0 in the velocities. Since O0y/0x, is itself a vector it follows that 
¢ij must be a (numerical) tensor, and it is immediately obvious that it can 
only be a multiple of 5;;. In fact if u, and v; are any two vectors, the associ- 
ated rational invariants are their squared lengths ua’, va” and the cosine of 
the angle between them, (wava)/[(ua?)(ve*)]”*. This fact shows that 
Caguag must be Ccbagua’g. Hence cy = cd, as was stated, and the force 
vector under consideration must be Oy/0x, up to a constant multiplier 
which may be absorbed into y. 
But the condition (5) upon this component of f; obviously yields 
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dy _ 
ds = 0, (Ai) 


i.e., the atomic potential remains constant along the world line of any particle 
of the perfect fluid. It is assumed that y vanishes along the free boundaries 
and in empty space. 

The formula for the corresponding atomic body force is, of course, 


PB. 
Sat = Se, (Ae) 


x4 


The primordially given atomic potential y supplies a useful mathematical 
instrument in the construction of a conceptual theory of matter and elec- 
tricity. Thus in 1927 (loc. cit.) I showed how the atomic potential might 
be used to obtain an atomic frequency equation which closely resembled the 
celebrated wave equation of Schrédinger. 

The model atom which I proposed was not positively unstable, although 
without rigidity. If it be required that the elementary constituents of 
matter, such as the proton, electron and neutron, cannot become locally 
concave this difficulty disappears (see my Oslo paper). We would have to 
suppose then that, in the moment when such concavity tends to be pro- 
duced, there arises a tensional normal force at the surface just sufficient to 
prevent it. Under such a condition a closely packed set of elementary 
constituents would necessarily have polyhedral forms, and this property 
would seem to indicate a possibility of crystalline structure and rigidity. 
There would obviously be a tendency of such everywhere convex bodies to 
maintain a roughly spherical form under collisions and other strong dis- 
turbances. 


V. THE ELECTROMAGNETIC POTENTIAL ¢, 


As indicated above, the electromagnetic force vector f,;, which arises from 
the terms linear in the velocities, is to be a vector of the form 


Since it is our intention to identify the potential y, with the usual electro- 
magnetic vector potential, we impose the condition 


re) (2e: ore re 
O%e Ne Ox) ote (2) 


where o is the density of electricity. Now the only possible constituent 
terms in fg, are to be obtained from 
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O%, 
— = %, 
OX, 


by a single contraction of indices and choice of a subscript i, in the three 
possible ways: 


ro Ova Ova 
— — — uy. 


ee ae ae 


Consequently the force vector in question is of the form 


0”; Ova Ova 
a Der, Ue +b vi Ug + ml Ut. 


But by the general requirement (5) this has to vanish fori = 1 when u = 1, 
Us = Ug = & = 0, Le., 


Se 


e+ He tte 0. 


Since there is no necessary relation between 0g;/0x; and 0y,/0*q, we infer 
thata +b=c=0. Hence this electromagnetic force vector is essentially 


ra) Oa 
far = 0 (32: _ a eg (Es) 


since the ponderomotive force is proportional to the electrical density co. 


VI. THE GRAVITATIONAL POTENTIAL hy 


By analogy with the Poisson equation in classical gravitational theory, 
we shall assume that a similar equation holds in each separate component 
of the symmetric energy tensor T;, and the corresponding component h,, of 
the symmetric gravitational potential h,,, namely, 


O*h4 
st = Sale G,) 





where the multiplier 87 is selected for reasons of convenience. This condi- 
tion (G,) is evidently the simplest analogous equation from the formal point 
of view. 

According to our initial hypothesis the types of terms which may enter 
in the corresponding gravitational force vector f¢;, arising from the quad- 
ratic terms in the velocities, are derived from 

Oh, 


— ut; 
Ox; 
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by a double contraction of indices and choice of the index 7. This yields 
six possible types: 








Ota Ohas hap 
ox, Ugg, Oa, Ugls, Oke U gui, 
Ohaa i aa 
ong U gui, One up", ae, upg’. 


But the last two types reduce respectively to Ohi_/O0%_ and Ohge/ Ox; since 
Ug” = 1, and so are of the degenerate type excluded by the hypothesis (c). 
Thus the most general available gravitational force vector is of the form 


Oia Oap 

aot << 

a om Ugg + 4 

But by (5) this must vanish fori = 1 if m4 = 1, m% = us = um = 0, so that 
Ohy Oat Ohaa 


rns te +e 


Since the quantities O/y/0%1, Ohe:1/O%_ and Ohee/O0x, are independent of 
one another we must havea + b =c=d=0. This yields for the gravi- 
tational force vector essentially the following expression 


ee 
Sot = p (2s Ox, Ug, (Ge) 


where the factor p is introduced since the gravitational force is proportional 
to the density p. 


Oha eee 
Uglsp +c = Ugu; + d axp Ugu; 


= 0. 











VII. THe COMPLETE THEORY 


We have then an energy tensor 7; given by (2) with p = p/2. The equa- 
tions of motion are given by (3) where 


fi — Sat + fri + fot 


The three terms on the right are defined by (Az), (E2), (Gz), and are of degrees 
0, 1, 2 in the velocities. Furthermore the equations (A1), (Zi), (Gi) deter- 
mine, respectively, the atomic, electromagnetic and gravitational poten- 
tials involved. Thus we have 21 dependent variables, p, o, us, ¥, 91, hay, 
and we have 20 equations in these variables, namely, the 4 equations of 
motion (3), the single equation (A;), the 4 equations (&,), the 
10 equations (G,) and in addition the relation u.2 = 1. This is as it 
should be since the electromagnetic potential ¢; is only determined up to the 
gradient of an arbitrary function, 0V/dx;. It is to be remembered that in 
empty space we take the atomic potential ¥ and the energy tensor T,, to 
vanish. 
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This system of equations is complete and consistent as a mathematical 
embodiment of matter, electricity and gravitation. The corresponding 
system in the generalized ‘theory of Einstein is incomplete to the extent 
that the equation of state of the homogeneous adiabatic fluid is not speci- 
fied, and it is inconsistent in that when two portions of the fluid collide, the 
equation of motion may break down. 


VIII. R&suMé OF THE NEW GRAVITATIONAL THEORY 


Let us now disengage as far as possible the new gravitational theory and 
consider to what extent its predictions are in agreement with the known 
facts. 

Matter is supposed to be either that special, mathematically satisfactory, 
homogeneous adiabatic fluid for which p = p/2 or, presumably, any form 
of matter in which the disturbance velocity is that of light under all circum- 
stances. 

We take 7;,; to designate the energy tepsor of matter and suppose that 
the equations of motion in the absence of a gravitational field (i.e., when 
only a small quantity of matter is present) may be written in the usual 
form: 


it... = 
Xa ies fo 





where fris a suitable force vector. It is then assumed that in the case of 
a gravitational field we may write the force vector f; in the form 


fi = fi + for 


where fg; designates the gravitational force vector.‘ 

We assume further that there is an associated symmetric gravitational 
tensor potential h,; such that a Poisson equation (G;) holds in each com- 
ponent of h,,; relative to T;;. In empty space 7; is taken to vanish. 

Under these circumstances fg; cannot be independent of the velocities 
since a proper force vector must be orthogonal to the velocity vector. 
Moreover since gravitational theory is reversible in time (in contradistinc- 
tion to electromagnetic theory), and since in the classical theory the gravi- 
tational force components are given by the components of the gradient of 
the gravitational potential, we assume that the gravitational force vector is 
homogeneous and quadratic in the velocity components, and homogeneous 
and linear in the first derivatives of the components of the gravitational 
potential, and furthermore that none of its components are degenerate, i.e., 
involve a factor u,2 = 1. 

In this way we obtain for the gravitational force vector fg, the unique 
expression (G:). Thus our complete system of equations is 
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OT 12 


orn 
OX i“ 5 = ST, tds 


OX.” 








fi t+ fow 


where we regard fg, and 7, as replaced by their explicit expressions. 

It is our intention in conclusion to indicate why this simple theory of 
gravitation set in the framework of flat (i.e., classical electromagnetic) 
space-time is in agreement with the observed facts. 


TX. GRAVITATION IN THE QUASI-STATIONARY STATE 


Let us suppose first that the portions of the perfect fluid are moving at 
small velocities relative to some frame of reference so that we have approxi- 
mately uw, = 1, uw = us = u% = 0. In this case we find in the coérdinates 
xs=t,2°=2,x* = y, x = 2, 


0? 0? 0? 

Bhy = —4xpby (a *.* dy? + 2) 
approximately so that hy, is negligible for 1 ¥ j while h,, fort = 1, 2, 3, 4 
reduce to the ordinary gravitational potential g. The gravitational force 
vector per unit of mass then reduces to the gradient of gz. Thus the theory 
is in first order agreement with the Newtonian theory. 


X. THE CENTRALLY SYMMETRIC STATE 


Suppose next that we have a sphere of the perfect fluid at rest with its 
center at (x, y,z) = (0,0,0). Ofcourse 7, and h,, will then be independent 
of the time ¢, and our Poisson tensor equation reduces to the equation writ- 
ten above in the same coérdinates so that we obtain for all 7 and j the exact 
equation outside of the sphere 


hy = - 543 (r, radial distance) 


where m is the mass of the fluid sphere. 
Thus the three exact equations of motion for a particle at (x, y, 2) at- 
tracted by the sphere are found to be of the type 


msx'r! 


2” 





x" = — ME BOE (alt + 9! + 5!8) + 


where the accent ’ indicates differentiation as to s. It is to be observed 
that the first term on the right yields the dominant Newtonian force com- 
ponents, the other two small terms being relativistic in origin. 

Now there is no essential restriction in assuming that the initial plane of 
the motion through (0, 0, 0) is the z-plane, whence we can at once conclude 
that z = Ofor all time. Hence we have only to solve the first two equations 
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with z = 2’ = 0,7 = /x? + y*. This is a readily integrable pair of equa- 
tions with x’, y’ and y, —x as two pairs of integrating factors. 
The differential equation of the path of the particle is seen to be 


22 (3) a ut) = gm + (gm +0) _ 7), 
dé 

where u = 1/r and 6 is the longitude; 4 and C are arbitrary constants of in- 

tegration. This may be integrated by an obvious quadrature. 

It is thus readily established that the resultant formulas for the advance 
of perihelion of the particle (x, y, 2) and the deviation of a ray of light 
(thought of as the path of a photon) in the gravitational field of the central 
sphere of matter are the same in their principal parts as in the Einstein 
theory. Furthermore the formula for the spectral shift toward the red is 
also in essential agreement with that theory. The exact expressions are, 
however, different.® 

Thus the simple theory of gravitation here outlined seems well adapted to 
explain the known physical facts. Like the Einstein theory, it has the ad- 
vantage of involving no arbitrary constants whatsoever. However, it is 
essentially different in that it presupposes. a framework of flat space-time 
instead of space-time curved by matter, and a basic form of matter in which 
the disturbance velocity is that of light. 


1“A Theory of Matter and Electricity,” ‘‘The Hydrogen Atom and the Balmer 
Formula,’’ these PROCEEDINGS, vol. 13, 1927. See alsomy article, ‘‘The Foundations of 
Quantum Mechanics,” in the Proceedings of the International Mathematical Congress at 
Oslo, 1 (1936). 

2 “Relativitatsprinzip und Gravitation,’ Physikalische Zeitschift, 13 (1912). 

3In a paper about to appear in the Revista de Ciencias de Lima, entitled “Sobre el 
Flidido Perfecto,” I have shown by direct integration that at least in two-dimensional 
space-time, such difficulties do not arise with the perfect fluid. 

‘It is a mistake to believe that the Einstein theory of gravitation does not similarly 
superimpose a gravitational force upon the other forces. It is only the mechanism of 
the superimposition which is different in the Einstein theory. In fact any physical 
theory without gravitation in flat space-time becomes one with gravitation in curved space- 
time when ordinary derivatives are replaced systematically by the covariant derivatives 
of the tensor calculus. 

5 The details of the new gravitational theory will appear in my article ‘El concepto 
de tiempo y la gravitaci6n” in the Proceedings of the Astrophysical Congress heid at 
Puebla and Tonantzintla, Mexico, in February, 1942. 
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GROUPS OF TRANSFORMATIONS OF THE NON-INVARIANT 
SUBGROUPS 


By G. A. MILLER 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ILLINOIS 


Communicated June 24, 1943 


All the non-invariant subgroups of a given group G are transformed under 
G according to a permutation group K whose degree is equal to the number 
of these subgroups, but not every permutation group of an arbitrary de- 
gree m is the group of transformations of the subgroups of some group. In 
particular, it will be proved in what follows that no imprimitive permuta- 
tion group is the group of the transformations of all the non-invariant sub- 
groups of some possible group and all the possible primitive permutation 
groups which are the groups of the transformations of the non-invariant 
subgroups of some group under this group will be determined. It therefore 
will remain only to determine some of the possible intransitive permutation 
groups which can be the groups of the transformations under a group of 
all its non-invariant subgroups of G. 

The operators of G which correspond to the identity of K are composed 
of all the operators of G which transform each of its subgroups into itself 
and they constitute a subgroup of G which has been called the norm of G.} 
To emphasize the fact that they constitute an important special character- 
istic subgroup of G we shall call this subgroup the character subgroup of G. 
A necessary and sufficient condition that the character subgroup of G is G 
itself is that G is either abelian or hamiltonian. In all other cases the 
character subgroup of G is a proper subgroup of G whose index under G is 
at least as large as 2, and when it is of this index it cannot be hamiltonian 
and it must therefore be abelian. 


To prove this fact it may first be noted that whenever the character 
subgroup of G is of prime index under G it involves all of the non-invari- 
ant subgroups of G. In particular, it involves all the operators of prime 
order contained in G. If the character subgroup of G is of index 2 under G 
it must therefore contain all the operators of order 2 which appear in G. 
There must be operators in G whose order is a power of 2 and which are 
not contained in the character subgroup of G. If this subgroup were non- 
abelian and would have an order which is a power of 2 it would also contain 
all the operators of order 4 which appear in G. If this were not the case an 
operator of order 4 which would not be found in the character subgroup of 
G would generate an invariant subgroup of G but its square could not be 
the square of the operators of order 4 found in the character subgroup of G 
since this subgroup contains all the operators of order 2 which appear in G. 

It therefore results that the operators of order 4 which would not appear 
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in the character subgroup of G would have different squares. As this is 
clearly impossible and since not all of the operators of G which do not ap- 
pear in its character subgroup could be of order 8 as a result of the proper- 
ties of the possible groups of order 16 it has been proved that whenever 
the character subgroup of G is of index 2 1t must be an abelian subgroup of G. 
When the character subgroup of G is of a larger prime index under G it is 
not necessarily abelian as results directly from the fact that a direct prod- 
uct of the quaternion group and a group which involves a prime number of 
conjugate subgroups, while all its other subgroups are invariant, contains 
a non-abelian character subgroup of prime index, where this prime is an 
arbitrary odd prime number. 

It is known that a regular group of any prime order is the group of the 
transformations of all the non-invariant subgroups of some group. We 
proceed to prove that no regular group of composite order can have this 
property. This will be done by showing that we arrive at a contradiction 
by assuming that a regular permutation group of composite order is the 
graup of the transformations of all the non-invariant subgroups of some 
group. Two facts which should be kept in mind in this proof are that all 
the non-invariant subgroups contained in.G would be conjugate under G 
and that all these non-invariant subgroups would appear among the opera- 
tors of G which would correspond to the identity of K. 

To simplify this proof we shall first consider the case when this regular 
group is of even order. To a permutation group of order 2 in K there 
would correspond in G a group whose character subgroup would be of 
index 2 under the group and hence this character subgroup would be abelian. 
Hence this group would have a commutator subgroup of order 2 and all of 
its operators which do not appear in the given character subgroup would 
generate this commutator subgroup since each of these operators would 
generate an invariant subgroup of G. If K would involve an operator of 
order 4 it may be assumed without loss of generality that this operator 
would generate the given permutation of order 2 contained in K. This 
clearly leads to a contradiction since the operators of G which correspond 
to the given permutation group of order 2 in K would transform all the 
non-invariant subgroups of G among themselves. This would also be 
true if K would involve a non-cyclic subgroup of order 4. 

The fact to be emphasized in this connection is that the given subgroup 
of order 2 in K would be invariant under K since the corresponding sub- 
group of G would be invariant under G. If the order of K is not divisible 
by 4 then K must contain a permutation of odd prime order which is 
commutative with the given permutation of order 2. To this permutation 
there would correspond in G an operator which would give rise to a com- 
mutator of odd order. As this is impossible we have arrived at a con- 
tradiction by assuming that the order of K is an even number. If K would 
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have an odd order it is easy to see that the character subgroup of G could 
not have an order which is a power of 2 and hence that this subgroup would 
again be abelian. It therefore results that K cannot be a regular group of 
composite order. 

If K would be an imprimitive non-regular group its subgroup composed 
of all its permutations which omit a given letter could not be maximal and 
hence G would contain non-invariant subgroups of different orders. These 
could not be conjugate and hence a proof of the fact that a non-regular 
imprimitive group cannot be the group of the transformations of all the non- 
invariant subgroups of some group under this group is a direct consequence of 
the fact that in such a group the permutations which omit a given letter 
form a non-invariant proper subgroup which is contained in a larger non- 
invariant proper subgroup of the imprimitive group. 

When K is primitive and non-regular its subgroups composed of all its 
permutations which omit a given letter omit only one letter and hence the 
number of these subgroups in K is the same as its degree. It therefore 
results that these subgroups are of prime order and hence that K involves 
a regular subgroup as an invariant subgroup. Moreover, this subgroup 
must be of prime order since the number of the non-invariant subgroups 
of K cannot exceed its degree. It therefore results that if a non-regular 
primitive group is the group of the transformations of all the non-tnvariant 
subgroups of some group under this group its order is the product of two dis- 
tinct prime numbers which are such that the larger of these two prime numbers 
diminished by one is divisible by the smaller. This condition is obviously 
both necessary and sufficient. 

The intransitive permutation group of prime order and of a degree which 
is an arbitrary multiple of this order is the group of the transformations 
of all the non-invariant subgroups contained in some group. In particular, 
it is clearly the group of the transformations under the group of all the 
non-invariant subgroups of the group obtained by forming the direct prod- 
uct of a prime power cyclic group and groups which satisfy the condition 
that each contains a number of conjugate subgroups which is equal to this 
prime number but contains no other non-invariant subgroup. It was noted 
above that such groups exist for every prime number and that when the 
prime is odd the group itself may be the group of the transformations of its 
conjugate subgroups, as is the case, for instance, when the group is the sym- 
metric group of order 6. 


1 Cf. Baer, R., Amer. Jour. Math., 61, 700 (1939). 
2 These PROCEEDINGS, 29, 106 (1943). 
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ON THE UNIFORM CONVERGENCE OF THE SOLUTIONS OF THE 
NAVIER-STOKES EQUATIONS 


‘By T. Y. THomas 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CALIFORNIA AT LoS ANGELES 
Communicated July 6, 1943 


Let ® be a fixed region of space completely filled by a viscous fluid. 
We assume § to be finite and to be bounded by one or more surfaces each of 
which may be at rest or in motion. -In the case of a fluid occupying an 
infinite region, for example, a hypothetical infinite pipe, §t, may be taken to 
be the periodicity cell which it is customary to introduce for convenience in 
such discussions. 

Denote by a#* (a = 1, 2, 3) and » the contravariant velocity compo- 
nents and the pressure function of a particular (stationary or non-stationary) 
fluid motion in ®. Similarly u* and will represent the corresponding 
functions for any other motion of the fluid in ®. It will be assumed that 
the viscosity u and density p of the fluid are constant and that all motions 
satisfy the Navier-Stokes equations and the equation of continuity. The 
quantities & = u* — #* and { = p — pare the components of an invariant 
€ called the disturbance. We shall derive general conditions under which 
all such disturbances converge uniformly to zero in §{ whereupon the above 
particular motion a p will be said to be stable relative to the class of mo- 
tions u* p (or disturbances £) in question. 


Consider the following inequality 


ap 8 (ES + Bye SF (ES + os, Pf’) = 0 (1) 


valid in any Riemann space (positive definite form). In (1) the quantities 
&¢, are the components of the covariant derivative of the spatial part of 
the disturbance ¢ and the functions f* are assumed to be continuous and 
have continuous first partial derivatives in the region ff. Expanding the 
left member of (1), integrating over Jt, making use of Green’s theorem and 
the boundary condition (* = 0 on the boundary of 9) we find 


S bap SES ERIV > S (fe — bff’) sap Pav. (2) 


In the case of periodic disturbances for which is the periodicity cell the 
condition of spatial periodicity combined with the above boundary condi- 
tion over the actual surfaces limiting § suffices for the derivation of this 
inequality. From (2) and the energy equation for the disturbance (de- 
rived from the Navier-Stokes equations) we obtain the inequality 
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T+ SWS — tesf"S?) tan + Dal PAVS0, — (3) 


in which the D.g are the components of the symmetric deformation tensor 
of the motion #* » and K is the kinetic energy of the disturbance. 


Now the functions f“ can be chosen so that the form 


(f5 — Sy8 FF gan & (4) 


is positive definite in the closure Jt of Rt, the condition for this being merely 
the existence of the condition ff — g,,f%f° > Oover R. Itis easily seen 
that this condition can be satisfied. For example if we introduce rectangu- 
lar Cartesian coérdinates x, y, z with origin in and take f! = F(x), f? = 
f*® = 0 the condition will be met if dF/dx — F* = c* where c is a constant 
different from zero. Integration of this equation gives F = c tan cx and 
hence F will be continuous and differentiable in $f if c is sufficiently small. 
It follows that (4) can be made positive definite in $i as above stated. 
Hence if the deformation tensor D is sufficiently small (at all times) the 
integrand in (3) will be positive definite. Under this condition it follows 
from (3) that K, the kinetic energy of the disturbance, will continually 
decrease and from the physical standpoint one must expect that this will 
result in the gradual “damping out” of the disturbance. A strict mathe- 
matical proof of the fact that the disturbance converges uniformly to zero 
can be made on the basis of certain assumptions concerning the con- 
tinuity and finiteness of the disturbance and certain of its derivatives and 
application of Ascoli’s theorem on the convergence of sequences. The de- 
tails of proof are essentially the same as those of a proof given in a previous 
paper: “Qualitative Analysis of the Flow of Fluids in Pipes,’’ Am. Jour. 
of Math., LXIV, 754-767 (1942). Understanding that these conditions 
are met we are thus led to the following result. Any motion ui p is stable 
in the region § relative to arbitrary finite disturbances provided that its 
deformation tensor D is sufficiently small. It follows in particular that if 
the deformation tensor vanishes the motion is always stable. 

As a specific example consider the stability of Poiseuille motion in a cir- 
cular pipe relative to arbitrary spatially periodic disturbances. For this 
motion the deformation tensor has the components Dis = Ds, = — wor/a? 
where we have referred the motion to cylindrical codrdinates r, ¢, 2 (the 
z-axis coinciding with the axis of the pipe) and where a is the radius of the 
pipe and wy the velocity of the fluid along the axis. The other components 
of the deformation tensor vanish. The conditions for the integrand in (8) 
to be positive definite can be expressed in the usual manner in terms of the 
elements of the symmetric matrix 
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ulf2 — ba ff?) 0 -" 
0 whe — @af?f) 0 
“> 0 ie tel T 


On observation we see that these conditions reduce to the single inequality 


fi — en Sf > %, (5) 


where R = woa/y is the Reynolds number, the quantity v being the kine- 
matical coefficient of viscosity. To treat the condition (5) let us make the 
substitution f, = —y,/y. The inequality (5) then becomes 


w+ Ty<0, W>0), 6) 


where the condition y > 0 is necessitated by the above requirement that the 
functions f, shall be continuous and differentiable in #. We can suppose 
the inequality sign in (6) to be replaced by the equality sign as a limiting 
case. Or if we like we can think of this equation as involving a slightly 
increased value of R in which case it will be equivalent for our purpose to 
the strict inequality (6). With this in mind we now suppose that y is de- 
pendent on r alone and make the change of independent variable x = r°/ 
whereupon the condition (6) becomes 


dy 
dr? 


A solution of (7) is given by the Bessel function Jo. In fact we have 


y= (qs) ~ (2) 


The requirement y > 0 in ® means that the argument of the above function 
J shall not exceed the value of the first root of this function in the interval 
Osr<a. This gives V4R/9 S 2.40r R S 13 (in round numbers). 
There seem to be two distinct mathematical problems concerning the 
above Reynolds number R according to experimental evidence. First, it is 
well known that large disturbances are damped out when R does not exceed 
a value around 1000. However, recent experiments have shown that even 
for very large values of R disturbances are likewise damped out provided 
they are sufficiently small.!_ There appears in fact to be no upper limit to 
the value of R beyond which small disturbances will not converge to zero. 
To account for this one may have recourse to the so-called linearized equa- 


+= “vt =v =0, Org. (7) 
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tions of the disturbance which can be justified strictly only under the as- 
sumption of arbitrarily small disturbances and this theory should thus lead 
to the stability of Poiseuille motion relative to such disturbances for any 
value of the Reynolds number. In spite of this, stability has been proved 
only for comparatively small values of R and even then only under the ad- 
ditional assumption of axial symmetry. Thus Orr’ has arrived at the value 
R < 180 for stability of Poiseuille motion relative to infinitely small axial 
symmetric disturbances. The problem of showing the existence of a criti- 
cal value of R at approximately 1000 is a problem involving finite distur- 
. bances and the general theorem of this paper has led only to the insuffi- 
cient value R S 13 for the stability of Poiseuille motion relative to these 
disturbances. As far as I am aware there is no previous discussion in the 
literature of the stability of Poiseuille motion relative to finite disturbances 
on the basis of the Navier-Stokes equations. 


1 See Prandtl, L., and Tietjens, O. G., Applied H'ydro- and Aeromechanics, McGraw- 
Hill, 32-35 (1934) for a more complete discussion and references to the literature. 

2 Orr, W. McF., “‘The Stability or Instability of the Steady Motions of a Liquid. 
Part II: A Viscous Liquid,” Proceedings of the Royal Irish Academy, (A), 27, 135 (1907). 

For a comprehensive discussion of stability relative to small displacements see Synge, 
J. L., “Hydrodynamical Stability,” Am. Math. Soc. Semicentennial Publications, I: 
Addresses, 227-269 (1938). 


RELATIONSHIP BETWEEN PARTICLE SIZE AND EFFICIENCY 
OF PALLADIUM-POLYVINYL ALCOHOL (Pd-PVA) CATALYSTS* 


By Louis D, Rampino,t Kevin E. KAVANAGH AND F. F. Norp 


DEPARTMENT OF ORGANIC CHEMISTRY,+ FORDHAM UNIVERSITY 


Communicated June 21, 1943 


Spectrophotometric measurements’ carried out on PVA and silver have 
proved that freezing is a suitable means for obtaining reproducible changes 
of the size of particles of colloids. Continued investigations in this labora- 
tory on colloidal noble metals, such as palladium, platinum and rhodium, 
in conjunction with synthetic high polymers and vanadium resulted in 
preparing supported catalysts, which have shown themselves to be ex- 
tremely efficient not only from a kinetic? but also from a chemical® stand- 
point. 

It was shown in a recent communication‘ that oriented nickel films with 
the same available surface as unoriented films had five times the activity 
as measured by ethylene hydrogenation. An attempt was made, ac- 
cordingly, to obtain insight into the causes for increased activity of Pd- 
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PVA catalysts by coérdinating ultracentrifugation of such solutions with 
measurements of hydrogen consumption, the rate of hydrogenation of 
nitrobenzene“ and determination of the size and shape of catalyst particles 
through electron micrographs. 

Measurements of Gas Consumption.—First, it was attempted to deter- 
mine whether any difference existed between the surface area or particle 
size of Pd in two types of colloidal supports. PVA was chosen as the repre- 
sentative member of the synthetic polymer series, and gum arabic (GA) 
as the representative of the naturally occurring type of colloid. The 
amount of gas adsorbed served as a means of measuring the surface area: 
for the surface area can be measured by the volume of gas adsorbed. Thus 
a difference in the volume of gas adsorbed by the Pd supported by PVA and 
by Pd supported by GA would serve as a means of distinguishing between 
the size of Pd particles in each medium. 

Apparatus.—A thermostat was set at 30°C. = 0.01. Special flasks were 
designed to fit constant pressure manometers, whereby the amounts of gas 
adsorbed were measured. The volume of the vessels which were cylindrical 
in shape was approximately 300 cc. Mercury served as the confining 
liquid in the manometers and dibutylphthalate was used in the supple- 
mentary capillary manometers. 

Two hundred and fifty cubic centimeter samples of 1% PVA solution 
containing 50 mg. of Pd were prepared from 125 cc. of a 2% PVA solution, 
120 cc. of distilled water and 5 cc. of a PdCl; solution containing 10 mg. of 
Pd per cubic centimeter, i.e., 1% with regard to Pd. Asa blank, a 1% 
PVA solution was made up from 125 ce. of 2% PVA and 125 ce. of distilled 
water. 

These solutions were immersed in the thermostat overnight to bring 
them to constant temperature. The next day, the special flasks, immersed 
in the thermostat and attached to the manometers, were evacuated by 
aspirator, filled with hydrogen and then flushed with Hz seven times. 

The hydrogen adsorption on the colloidal palladium was the difference 
between the PVA-Pd solution and the blank. The blank reading was 
added if the volume in the blank solution increased ; subtracted if the volume 
in the blank solution decreased. 

The total. gas consumption figure thus obtained was diminished by the 
amount of hydrogen which is necessary to reduce the PdCl, (20 mg. Pd) 
to metallic Pd. Thus, the absolute adsorption is the total gas consump- 
tion figure minus 4.20 cc. For comparison purposes, however, the total 
gas consumption figure was used. 

The results obtained with PVA-Pd and GA-Pd, in which the colloids 
were present in a concentration of 1% and 0.5%, respectively, are pre- 
sented in table 1. 
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TABLE 1 


HypRoGEN CoNsuMPTION BY GA-Pd anp PVA-Pd Soxutions (100 Cc.) WuicnH Have 
Stoop In THERMOSTAT 16-20 Hours 


cc. oF Hz ADSORBED 


COLLOIDAL CATALYST 0.5% 1% 
GA-Pd 8.30 + 0.20 8.30 = 0.20 
PVA-Pd 6.50 = 0.20 3.90 = 0.20 


Two observations are of interest. First, a greater hydrogen consumption 
indicating a greater surface, is possessed by the inferior GA-Pd catalyst; 
further, an increase in consumption is shown by the PVA series when the 
concentration is decreased to 0.5%, whereas the GA series shows identical 
hydrogen consumption for both concentrations. 

It had been noticed that the reddish colloidal solution containing PVA 
and PdCl, had turned to a blackish solution on standing for several hours 
(without hydrogen) in the thermostat. The GA-PdCl, solution did not 
exhibit this property but retained its light yellow translucent color. 
To circumvent this darkening, the PVA solutions were made up, let stand 
in the thermostat overnight and then on the next morning, immediately 
before introduction into the measuring vessels, made up to volume of 
250 cc. by the addition of 5 cc. of the PdCl, solution. The same was 
done with the GA catalyst. By this treatment, no darkening occurred in 
the PVA solutions. 

With the GA-Pd solutions, no change in hydrogen adsorption was noted. 
However, striking results were obtained with the PVA catalysts. PVA 
(1%)-Pd showed a total hydrogen consumption of 8.05 cc. + 0.20; while 
a PVA (0.5%)-Pd showed a gas consumption of 8.25 cc. = 0.20. The 
values given here as well as all other figures previously recorded, were the 
mean values of six measurements. It was apparent, then, that there was 
no significant difference between the adsorption of 1% and 0.5% PVA 
when the solutions were made up immediately before placement in the 
measuring vessels. These results are recorded in table 2. In table 3, all 
values are converted from ‘“Total Gas Consumption”’ to “Absolute Adsorp- 
tion,” i.e., the amount of hydrogen required for reduction of PdCl, is sub- 
tracted. 


TABLE 2 


H, CONSUMPTION BY ‘‘FRESHLY PREPARED’”® PVA-Pd anp GA-Pd Catatysts. 100 
Cc. or SOLUTION CONTAINING 20 Mc. or Pd ARE USED IN ALL CASES 


cc. oF Hz ADSORBED 


CATALYST SOLUTION 0.5% 1% 
GA-Pd 8.30 = 0:20 8.30 += 0.20 
PVA-Pd 8.25 = 0.20 8.05 + 0.20 


* “Freshly prepared” is used to designate solutions which had the PdCl, added im- 
mediately prior to the placement in measuring vessels. 
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TABLE 3 


ABSOLUTE Hz ADSORPTION BY “STANDING’” AND “‘FRESHLY PREPARED” PVA-Pd AND 
_ GA-Pd Catatysts 


cc. H: ADSORBED 


CATALYST SOLUTION 0.5% 1% 
GA-Pd (standing) 4.10 = 0.20 4.10 = 0.20 
GA-Pd (fresh) 4.10 = 0.20 4.10 = 0.20 
PVA-Pd (standing) 2.30 = 0.20 0.00 
PVA-Pd (fresh) 3.85 + 0.20 4.05 = 0.20 


° “Standing” is used to designate solutions wherein the PVA and PdCl, stood in con- 
tact at least 17 hours before transfer was made to the measuring vessels. 


On analysis of these figures it is evident that the adsorptions observed 
in PVA-Pd solutions which are “freshly prepared,” approximate those 
shown by GA-Pd. However, in ‘standing’? PVA-Pd solutions a much 
smaller hydrogen adsorption than in GA-Pd was measured, and the 1% 
“standing” PVA-Pd showed no actual hydrogen adsorption at all, when 
allowance was made for the amount of hydrogen required for reduction. 

The explanation of these anomalous results was that the PVA upon 
standing in contact with PdCl, caused.a reduction of the PdCl, to 
metallic Pd. This would make intelligible the acquisition by the col- 
loidal solutions on standing of a blackish shade. It appears from the data 
that 1% PVA is able to accomplish this stoichiometrically whereas 0.5% 
PVA does not fully reduce PdCh, only effecting a partial reduction in the 
overnight time interval. 

This explained the apparent non-adsorption of hydrogen by ‘‘standing”’ 
PVA (1%)-Pd solutions, while a less active catalyst GA-Pd exhibited an 
adsorption of 4 cc. 

Thus there was no significant difference in hydrogen adsorption displayed 
by PVA-Pd or GA-Pd solutions. Therefore, by this procedure no differ- 
ences in particle size or available surface could be established between 
GA-Pd and PVA-Pd to explain the superiority of the PVA-Pd catalyst. 

Additional proof for the fact that PVA reduced the PdCl, to Pd is af- 
forded by the observation that upon determining the pH of the PVA- 
PdCl, solution immediately after mixing and then after standing, an in- 
crease in acidity was noted, due to the formation of HCl by the reduction 
of PdCh. 

Ultracentrifugation and Rate of Hydrogenation.—Meanwhile, the prob- 
lem of the relationship between particle size and rate of catalytic hydrogen- 
ation had been attacked from a different angle. The method of approach 
in this instance was to ultracentrifuge colloidal PVA-Pd solutions at various 
speeds, thus effecting a separation of the colloid into small Pd particles 
which would remain in the supernatant liquid and into the large particles 
which would be precipitated by the centrifuging. By quantitative deter- 
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mination of the Pd present in supernatant liquid and in residue and then 
by measuring the hydrogenation rates of equivalent quantities of Pd from 
both, an insight could be obtained into the effect of particle size and avail- 
able surface on the activity of the catalyst. Through the courtesy of Dr. 
R. W. G. Wyckoff, of Lederle Laboratories at Pearl River, New York, 
PVA (1%)-Pd solutions were ultracentrifuged at 200 and 500r. p.s. There 
was a noticeable difference in the shade of the two supernatant liquids, 
that which was centrifuged at the higher speed being much lighter. 

The supernatant liquids were analyzed in the following way: A 1% 
PVA solution, which had been centrifuged at 500 r. p. s., was analyzed to 
determine the quantity of PVA in the top one-third and in the bottom one- 
third. Two cubic centimeters of the top one-third contained 18.02 mg. 
of total solid. The bottom one-third contained 20.50 mg. of solid per 2 cc. 
of solution. This gave a total solid in the two portions of 38.52 mg. per 
4 cc. instead of the theoretical 40.00 mg. for a 1% solution. The total solids 
were then ignited for an hour to give 1.76% ash. This percentage of ash 
was then used to determine the quantity of Pd present in 5 cc. of super- 
natant liquid from the PVA-Pd solutions centrifuged at 200 and 500 
r.p.s. Five cubic centimeters of supernatnat liquid were taken to dryness 
and weighed to give total solids. Total solids were ignited and weighed to 
give ‘‘ash.”’ This ‘‘ash’’ was in fact ash plus Pd. The weight of ash present 
could be obtained by multiplying the total solids by the percentage of ash 
previously found, i.e., 1.76%. The weight of Pd could then be determined 
by substracting the weight of the ash thus obtained from the weight of Pd 
plus ash. 


TABLE 4 


DETERMINATION OF MG. oF Pd PRESENT IN 5 Cc. OF SUPERNATANT LIQUID AT DIFFER- 
ENT RATES OF ULTRACENTRIFUGATION 


ASH IN MG. 
(TOTAL Pd In mo. 
TOTAL Pd + asa % ASH SOLID X (asH + 
SAMPLE SOLID IN MG. IN PVA % Asx) Pd)—asu 
5 cc. (S. S.)* at 
200 r. p. s. 44.43 1.73 1.76 0.78 0.95 
5 ce. (S. S.) at 
500 r. p. s. 43.94 1.28 1.76 0.77 0.51 


* S. S. stands for supernatant solution. 


From the data it is seen that there was twice as much Pd present in the 
supernatnant liquid of the sample ultracentrifuged at 200 r. p. s. as there 
was in the sample centrifuged at 500 r: p.s. Using these data, it was pos- 
sible to set up a kinetic experiment in which the hydrogenation rate of 
equivalent quantities of Pd of different particle size could be compared; 
it is self-evident, of course, that the size of the Pd particles in the sample 
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centrifuged at 500 r. p. s. was smaller than that of the particles in the 
sample centrifuged at 200 r. p.s. If volumes of the supernatant liquids of 
the two centrifuged samples were taken which contained equivalent quanti- 
ties of Pd, the hydrogenation rate of the supernatant liquid from the sample 
centrifuged at 500 r. p. s. should be greater than that of the sample cen- 
trifuged at 200 r. p. s. because of the greater surface of the former. This 
assumption, which correlates catalyst efficiency and particle size, proved 
correct. 

The hydrogenation rates obtained with these two catalyst solutions con- 
taining equivalent quantities of Pd, ie., 2 mg., are shown in figure l, 
where it is demonstrated that the larger surface area in the 500 r. p. s. 
centrifuged sample is responsible for a greater efficiency of the catalyst. 

Next the residues from the ultracentrifuging were analyzed to determine 
their Pd content so that aliquots containing 2 mg. of Pd could be taken and 
then compared with the activity of the particles of much smaller size in the 
supernatant liquids whose hydrogenation rates were just established and 
plotted in figure 1. 

The results of analyses are recorded in table 5. 


TABLE 5 
32-cc. SAMPLE Pd IN VOL. OF MG. Pd PER TOTAL MG. MG. Pd IN 
CENTRIFUGED 32 cc. S. S.e 5 ce. S. S. Pd in S. S. RESIDUE 
200 r. p.s. 12.48 27.0 1 5.40 7.08 
500 r. p. s. 12.48 26.4 0.5 2.64 9.84 


2 “S_ S.”’ is used to designate the supernatant solution. 


In order to conveniently apply the sediments, 5 cc. of distilled water were 
added to each bringing the volumes up to 6.3 cc. and 5.1 cc. for the 500 
r. p. s. and 200r. p.s., respectively. Thereby, 9.84 mg. of Pd and 7.08 mg. 
of Pd were present in 6.3 cc. and 5.1 cc., respectively. By calculation, 
the proportional volumes containing 2 mg. of Pd were determined and it 
was found that 1.28 cc. of the 500 r. p. s. sediment solution and 1.44 cc. 
of the 200 r. p. s. sediment solution were required. These volumes were 
then made up to 40 cc. in the same manner as were the supernatant liquids. 

From figure 1, a convincing view of the effect of particle size on the hy- 
drogenation rate may be obtained. For therein is presented the hydro- 
genation velocity toward a common substrate of analogous colloidal solu- 
tions of PVA-Pd catalysts containing equivalent quantities of Pd, viz., 2 
mg. The difference between the rates of supernatant liquid and sedi- 
mented residues is striking. In the 500 r. p. s. sample, the sedimented ma- 
terial is about 13 times less active than the supernatant liquid using the 
adsorption of hydrogen at the end of five minutes as a basis for comparison: 
using fifteen minutes’ adsorption as a basis, the sedimented material is 
about 7'/: times less active. With the 200 r. p.s. sample, the sediment 











252 CHEMISTRY: RAMPINO, ET AL. Proc. N. A. S. 


is about 5°times less active than the supernatant using the adsorption at 
the end of five minutes as a basis of comparison; using fifteen minutes’ ad- 
sorption as a basis, it is about 4!/; times less active. Certainly, these 
ratios may only be attributed to the differences in surface area of the par- 
ticles. 

It was understood that the concentrations of PVA in the residue and in 
the supernatant solution were different, viz., 25 mg. and 90 mg., respec- 
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FIGURE 1 
Comparison of Supernatant Liquids and Residues 

Substrate: 2 mg. Pd in 40 cc. of 50% alcohol. Acceptor: 0.5 cc. 
C.,H;NO;. A, supernatant liquid of sample centrifuged at 500 r. p. s. 
B, supernatant liquid of sample centrifuged at 200 r. p.s. C, residue 
of catalyst centrifuged at 200 r. p. s. D, residue of catalyst centri- 
fuged at 500 r. p. s. All measurements are converted to 25°C. and 
760 mm. 


tively, in the 200 r. p. s. centrifuged sample and 21 mg. and 175 mg., 
respectively, in the 500 r. p. s. centrifuged sample. From the data ob- 
tained by Rampino,' it can be seen that the effect of PVA concentration 
would be insufficient to explain the striking differences in the measured 
efficiencies of the palladium in the residue- and supernatant-catalyst. 
The superiority of the supernatant liquid of the 500 r. p. s. centrifuged 
sample to the supernatant liquid of the 200 r. p. s. sample is also noteworthy. 











FIGURE 2 
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FIGURE 5 
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It appears from these results, that the bulk of catalytic activity rests pre- 
ponderantly with the smaller particles. Thus surface area and particle size 
play a paramount réle in hydrogenations performed with PVA-Pd catalysts. 

Correlation of Electron Micrograph’ Measurements and Kinetic Experi- 
ments.—Since no significant differences in particle size could be established 
by the gas consumption experiments previously described, it was deemed 
necessary to resort to a more accurate and direct means of measuring the 
particle size, viz., by applying the electron microscope. 

Figure 2 represents an electron micrograph of PVA alone. This showed 
a tendency to crystallize in acicular crystals of from 0.03 to 0.4 micron in 
length; the shape and density of these crystals was very different from 
that of the Pd seen in later micrographs. This micrograph was obtained 
from 0.5% PVA solution magnified 30,000 times. 

In figures 3 and 4 two PVA-Pd catalysts are shown. They are identical 
in composition: they are both 0.5% in PVA (Du Pont’s-Rh-391) and con- 
tain 10 mg. of Pd per 50 cc. of HxO. They only vary in the method used 
to reduce the palladium chloride to metallic Pd. The catalyst shown in 
figure 3 has been reduced by hydrogen, i.e., by introducing the aqueous 
PVA-PdCl, solution into a vessel and shaking for five minutes. The cata- 
lyst shown in figure 4 has been reduced by PVA itself, i.e., by letting the 
PdCl, stand in contact with the PVA whereby the slow reduction of the 
PdCl, to Pd, takes place. This micrograph, with its characteristic Pd par- 
ticles clearly visible, conclusively proves the fact that PVA reduces PdCl, 
to Pd, which was proposed to explain the anomalous gas adsorption results. 

In these two micrographs, the black dense metallic Pd particles are 
clearly visible in hexagonal crystals, in which allotropic form palladium is 
known to exist.6 By far the greater number of particles assume an ultra- 
microscopic hexagonal configuration but some triangular shaped crystals 
can be noted. In these micrographs, the magnification is also 30,000 di- 
ameters. In both of these, there is no evidence of the PVA with its char- 
acteristic acicular crystals. 

The possibility had been considered that these two PVA-Pd catalysts 
might have different hydrogenation efficiencies. The different modes of 
reduction could have produced Pd particles of varying sizes and of different 
surface areas and consequently two catalysts of varying effectiveness 
could have been formed. This seemed probable since one method of reduc- 
tion, viz., by means of hydrogen, was of short duration, about five minutes, 
while the other, viz., by means of PVA, took in the case of 0.5% PVA, 
more than twenty hours. Kinetic experiments were run in an attempt to 
establish whether this difference existed. However, no significant difference 
could be established. Both catalysts gave hydrogenation rates of the same 
magnitude, within the experimental error, using nitrobenzene as an ac- 
ceptor. 
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A statistical evaluation, wherein a size frequency curve of the single 
particles is determined from the electron micrographs of these two catalysts 
provides an explanation of their identity in kinetic behavior. A size fre- 
quency curve on the micrograph in figure 3, i.e., the PVA-Pd catalyst in 
which the PdCl, was reduced to Pd by hydrogen, gave a ds size” (average 
particle size with respect to surface) of 0.07 micron. The same type of 
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FIGURE 6 


Comparison of “Supernatant” and ‘“‘Normal’’ Catalysts 
Substrate: 2 mg. Pd per 40 cc. of 50% alcohol. Acceptor: 0.5 cc. 
C.H;NO:. A, supernatant liquid of PVA-Pd centrifuged at 200 r. p.s. 
B, normal H2-reduced PVA-Pd. All measurements are“converted to 
25°C. and 760 mm. 


analysis on the figure 4 micrograph, i.e., the PVA-Pd catalyst in which 
the PdCl, was reduced to Pd by PVA itself, gave a ds size of 0.08 micron. 
There is only a slight difference in particle size of these two differently re- 
duced catalysts, and this difference could not be expected to reveal signifi- 
cant deviations in the hydrogenation velocities produced by the two cata- 
lysts. Once again a relationship between surface area and catalyst effi- 
ciency is established. 

Figure 5 is an electron micrograph, at a magnification of 30,000, of the 
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supernatant liquid of the PVA-Pd catalyst which was ultracentrifuged at 
200 r. p.s. This micrograph was not made until a year and a half after the 
sample had been prepared and ultracentrifuged. As would be expected of 
a colloid which had been standing for such a length of time, a great deal of 
clustering had taken place, which is easily discernible in the micrograph. 

However, the presence of very fine grain Pd is unmistakable. An esti- 
mate of the ds size amounts to 0.009 micron. This value might have to 
be revised slightly upward because the larger particles cannot be ac- 
curately measured in the clusters. There can be no doubt, moreover, that 
the particle size of much of the Pd is smaller than in the previous micro- 
graphs. Using the approximation of 0.009 as the dy size, the particle size 
seems to be about one-eighth that of the PVA-Pd Hz: reduced catalyst, 
wherein the ds size was shown to be 0.07 micron by a size frequency curve. 

It has been previously described how a kinetic experiment was run on the 
supernatant solution of this 200 r. p. s. sample. For comparison purposes, 
the hydrogenation velocity of a non-centrifuged H,-reduced PVA-Pd cata- 
lyst, containing 2 mg. of Pd and of similar composition to the catalyst pre- 
pared from the supernatant solution of the 200 r. p. s. sample, was meas- 
ured. 

These two hydrogenation velocities are plotted in figure 6. 

It can be seen that the ‘‘supernatant”’ catalyst containing the same quan- 
tity of Pd as the analogous ‘‘normal” catalyst is more than twice as ef- 
ficient. The analogous ‘normal’ catalyst had a particle size about 8 
times larger than that of the ‘‘supernatant’’ catalyst; 0.07 micron (statis- 
tically determined) to 0.009 micron (estimated). From the results of these 
experiments, it is evident that there is a definite relationship betwee1. “he 
surface area of the catalyst and rate of hydrogenation. 

The only two organic colloidal catalysts to be investigated by von 
Ardenne® were Pt-Paal (sodium lysalbinate) and Pt-Skiia (GA). He esti- 
mates from his micrographs that the particle size in the former is about 
0.003 to 0.01 micron and in the latter, it averages about 0.025 micron. Of 
course, these measurements are made on platinum catalysts but it is rea- 
sonable to assume that the Pd particle size in these colloids would be of a com- 
parable magnitude. Ifsuch an assumption is valid, certain conclusions may 
be drawn. The particle size of Pt and Pd, too (by analogy), in these two 
colloidal catalysts may be smaller, according to von Ardenne’s measurements, 
than the particle size in PVA-Pd as determined from our micrographs; 
the Pt-Skita (GA) having a particle size diameter of 0.025 micron (esti- 
mated) to the 0.069 micron (calculated) diameter of PVA-Pd particles. In 
spite of this, the superiority of the PVA-Pd catalyst has been established. 
Also to be recalled is that no difference in particle size could be revealed by 
the gas adsorption measurements previously described. Consequently, 
although the preponderant effect of the particle size on the rate of hydro- 
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genation is definitely established, there seems to be another factor, a 
quality factor, the extent of which is, so far, not estimable, which is in- 
dependent of particle size and which makes the synthetic high polymer 
carrier far superier to other colloidal carriers, e.g., GA, even though the 
particle size in the two types of carriers may be the same. 

Whether the striking formation of ultramicroscopic hexagonal and tri- 
angular shaped Pd metal crystals in the presence of PVA has bearing on 
the original cause of the effectiveness of this type of catalyst, could be, 
perhaps, decided by comparing electron micrographs of catalysts deposited 
on solid bodies. 

Summary.—By coordination of measurements of hydrogen adsorption, 
rates of hydrogenation and determination of the size and shape of catalyst 
particles it was established that the preponderance of catalytic activity 
lays with the Pd particles of smaller size. However, the superiority of the 
synthetic polymer supported catalyst over the GA type is to be attributed 
to an additional quality factor. 


* This work gained the support of grants of the Bache Fund of the National Acad- 
emy of Sciences, and the Penrose Fund of the American Philosophical Society, and is 
abridged from a part of the dissertation of K. E. K., submitted in partial fulfilment of 
the requirements for the degree of Doctor of Philosophy to the Graduate Faculty of 
Fordham University. 

¢ Present address: Shell Development Co., Emeryville, Calif. 

¢ Communication No. 30. 

§ Our appreciation is expressed to Mr. Ernest F. Fullam of Interchemical Corporation 
of New York for his courtesy in preparing numerous electron micrographs and a discus- 
sion of the results, and to Baker and Co., Inc., Newark, N. J., for donating the noble 
metal salts used in this investigation. 


1 Leichter, H., and Nord, F. F., Biochem. Z., 295, 226 (1938). 

2? Rampino, Louis D., and Nord, F. F., Jour. Amer. Chem. Soc., 63, 2745, 3268 (1941). 

8 Kavanagh, Kevin E., Ibid., 64, 2721 (1942); Rampino, Louis D., and Nord, F. F., 
Ibid., 65, 429 (1943). 

* Beeck, O., Smith, A. E., and Wheeler, A., Proc. Roy. Soc. (London), A117, 62 (1940). 

4° Nord, F. F., Ber., 52, 1705 (1919). 

5 Dissertation, Fordham University, 1941, figure 5. 

® Naumann, C. F., and Zirkel, F., Elemente der Mineralogie, 15th edition, Leipzig, 
Engelmann, 1907, p. 421. 

7 Green, H., Jour. Franklin Inst., 204, 713 (1927). 

8 yon Ardenne, M., and Beischer, D., Angew. Chem., 53, 103 (1940). 








VoL. 29, 1943 BACTERIOLOGY: WILSON AND WORCESTER 257 


THE DETERMINATION OF L.D.50 AND ITS SAMPLING ERROR 
IN BIO-ASSAY, III 


By EpwIn B. WILSON AND JANE WORCESTER 
HARVARD SCHOOL OF PUBLIC HEALTH 


Communicated June 24, 1943 


1. The discussion of this problem has been conducted on the hypothesis 
that the number 1 of animals used in each of the dilutions is the same. This 
restriction is convenient in that it leads to the possibility of preparing a 
table for the determination of the constants and of their standard errors 
when there are three dilutions.1_ Formulae can, however, be developed for 
the more general case where the numbers are not the same, but ™, me, m3. 
The rule of maximum likelihood leads to the equations: 


L 
< = —ad(2s, — m) + adn, tanh a(x; — y) = 0, (1) 
uf 
oL 
in = (2s, — m)(x, — ¥) — U(x, — 7) tanh a(x, — y) = 0. (2) 


Placing x1, %2, x3 equal to x — c,x,x + ¢ and setting the derivatives equal 
to zero, introducing X = tanh a(x — y) and C = tanh ac we have the 
simultaneous equations for X and C 


C2[m2X* — AX? + (m + 3)X] + C(ns — m)(X? —1) +A — NX =0, 
(3) 


C?[(n3 - m)X —_ 2BX?] + C(ns a m)(X? = 1) -}- 2B — (ns —_ m1) ye 0, 
(4) 


where 


A = 2(s; + Se + 53) — N, 2B = 2s3 — 25, — m3 + m (5) 


and N = m, + m + nz is the total number of animals used. Here A and B 
which take the place of A and B in the earlier treatment are in terms of 
numbers, not of proportions, of the animals affected. The elimination of C 
gives a cubic equation for X only slightly more complicated than that 
found before, viz., 


Nn,X* — [A(N + m) — 2B(n3 — m)]X? 
+ [A? — 4B? + N(ms + m) — (ms — m)?7)X 
im: [A(ns +4 nm) i 2B(n3 — m) | an QO, 
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Having solved this numerical equation for X, one may obtain C from the 
equation 


Ae —A(n3 — m) + 2B(ns + m) + mo(m3 — m) X 
No(N3 oe m)X? —_ [A(ns3 _ n}) — 2B(n3 —_ m) |X aa Anns. 
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and then 


a= 


: 
c 


where the functions tanh~' may be replaced by their equivalents in 
logarithms. 


1 
tanh-! C, y =x — —tanh-'! X, 
: Qa 
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The sampling errors may be obtained from 





OL 
— = —a*In, sech? a(x, — y) = —4a72n,P,Q,, 
oy? 
OL 
3a? => — In x, —_ 7)? sech? a(x, ~— 7) = —42nP Ox, —_ 7)’, 
OL 

= aln,(x, — y) sech? a(x, — y) = 402 P,Qi(x; — 1), 
Oady 

when P, are the fitted values, by 
a 1 OL e 1 OL 


“ne. eae 
where H is the Hessian of L with respect to y and a. Asa matter of fact 
these formulae are valid for any number of dilutions provided a solution 


for a and y and consequently for P; upon the curve has been found in any 
manner, presumably by successive approximation from the equations 


a] 2), |, OL EL, Om 
Ovo «Oy? dad7\o wala Valo Dady da? 





by + by + ba = 0. 




















0 0 0 


2. In case one is working with a sufficient number of dilutions and in 
case the observed values P; do not fluctuate too much from a straight 
line on growth paper one may get an approximate solution by graphical 
methods. For example, if we have the following logarithmic doses x and 
proportions of population affected? 


x 1.755 1.784 1.811 1.837 1.861 1.884 
n 28 27 30 31 30 29 
s 9 14 23 29 29 29 
> i 0.329 0.519 0.767 0.936 0.967 1.000 


we may plot the first 5 points (P' = 0 or 1.0 cannot be plotted) and draw 
the line. The 50% end-point is read as at 1.777. The fitted growth curve 
being 

tanh—!(2P—1) 


= !1 1 — = 
P; /s + /s tanh a (x 1.777), a x — 1.777 





permits the determination of a from any values P; and x read from the 
line. At x = 1.837, P; = 0.935. Hence a = 22.2 and P; = 1/2, + 
1/, tanh 22.2 (x — 1.777) is the line as determined graphically. 

3. The values of P; as fitted may be read* from the line and the calcu- 
lation form may be filled out. 
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x 1.755 1.784 1.811 1.837 1.861 1.884 
P, 0.280 0.587 0.810 0.9385 0.975 0.992 
Qy 0.720 0.413 0.190 0.065 0.025 0.008 
P; —P —0.049 0.068 0.043 —0.001 0.008  —0.008 
ys —1.37 1.8 1.29 -0.08 0.2%  -—0.23 
x-—y¥ —0.022 0.007 0.034 0.060 0.084 0.107 
(s — s)(x — vy) 0.0301 0.0129 0.0489 -—0.0018 0.0202 —0.0246 
4nP,Q, 22.6 26.1 18.5 7.5 2.9 0.9 
4nP,Q,(x — +) -—0.50 0.18 0.63 0.45 0.25 0.10 
4nP,Q,(x — y)? 0.011 0.001 0.021 0.027 0.021 0.011 
Now, ‘ 
oL 
—| = 2aX(s; — s) = 77.8, ~ = —22(s, — s)(x — vy) = —0.161, 
Ov|o alo 
2 2 2 
oF = —38,700, = = —0.092, = = 24.6. 
O77! 0a*\o ady 





Then 77.3 — 38,7005y + 24.66a = 0, —0.161 + 24.65y —0.0925a = 0 
give as solutions dy = 0.00106, a = —1.47 and y = 1.77806, a = 20.73 
as a first arithmetic approximation subsequent to the graphical fit. If 
the graphical fit be considered an adequate approximation for computing 
o, and o, we have for the Hessian H = 2955 and hence 

38,700 0.092 


2 = —— = 131, 2 = — = 0.000031, 
¥ 2955 sie Sige 


which gives ¢, = 3.6, «, = 0.0056 and we write‘ 
vy = 1.7781 = 0.0056, a = 20:7 = 3é. 


In this example the logarithmic doses x did not advance by uniform 
amounts so that it is impossible satisfactorily to determine a and y from 
the table which was computed on the basis of an equal advance in x and 
constant values of m; but if one should ignore these points and take the 
lst, 3rd and 5th observed P’s we should find A = 1.126, B = 0.638 giving 
“yy” = —0.63 and ‘‘a” = 1.00. With c = 0.053 we have y = —0.033 and 
a= 19. The 50% end-point would be at 1.811 — 0.033 = 1.778. The 
values 1.778 and 19 in place of 1.777 and 22.2 determined graphically or 
1.7781 and 20.7 obtained by the first arithmetic approximation are much 
the same thing in view of the standard errors estimated from the graphical 
solution. We could not expect to get the values of the standard deviations 
from the table because they would be estimated from only about one-half 
of the observations and furthermore the estimates depend on the actual 
points used on the growth curve as well as on the number of observations; 
in fact a calculation gives 1.778 + 0.0080 and 19 + 4.2 in place of 1.7781 


Il 
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+ ().0056 and 20.7 + 3.6. In this case, however, the use of the table and 
only three of the six observed proportions gives a fairly good answer. 

4. There need be no restriction to three points or to identical values of 
n, If in (1) and (2) we write x, = x + ic, assuming the logarithmic doses 
to proceed in a uniform manner (or the dilutions to be a geometric progres- 
sion) we have for the determination of y and a 


—aX(2s, — m) + adn, tanh a(x + ic — y) = 0, 
(2s, — m) («@ + ic — y) — Uni(x + ic — y) tanh a (x + ic — y) = 0. 


If the first equation be multiplied by (« — y)/a and added to the second, 
the latter becomes 


r1(2s; a n:) = Lin, tanh a (x “+ 1c — 7) = Q. 


The equations will be algebraic in X = tanh a(x — y) and C = tanh c but 
of too high an order to be readily used. It will, however, be noticed that 
the solution depends on the particular two constants 


A = X(2s, on ni), 2B = ~~ (2s, ai nN). 


For this reason, if some standard set-up involving more than three dilu- 
tions in geometric progression were to be adopted for wide-spread use it 
would be possible, and might be well, to compute a double entry table for® 
y and a and their standard deviations in terms of A and B. 


1 These PROCEEDINGS, 29, 79-85, 114-120, 150-154, 207-212 (1943). 

2 For growth paper, see Wilson, E. B., these PROCEEDINGS, 11, 451-456 (1925). The 
data used in the illustration are taken from Bliss, C. I., p. 154 of Ann. Applied Biol., 22, 
134-167 (1935) or from Irwin, J. O., p. 17 of Suppl. Jour. Roy. Stat. Soc., 4, 1-60 
(1937). Those authors fit only the last six of the eight points. In the figure the limits 
of P at one standard deviation as determined from the formula on p. 209 of Jour. Amer. 
Statis. Assoc., 22, 209-212 (1927) have been inserted and two lines have been drawn, one 
fitting (graphically) the last six and the other all eight points. It is clear that the 
graphical fit to the eight points is not bad; indeed if one actually fits the eight points 
arithmetically one finds y = 1.7714 + 0.0056, a = 17.1 + 2.1 with a x? value of 3.35 
which for six degrees of freedom gives P = 0.76 with a probability that because of the 
small numbers in some of the cells, the true value of P would be even larger. It is true 
that the deviations seem to be of such signs as to indicate curvature of the observed 
points to the fitted line, but in view of the small value of x? and high value of P one 
might hesitate to draw that inference and to reject the first two points from this statis- 
tical evidence. The real reason for the rejection of these two points lies in the fact that 
in other corresponding sets of data the mortalities at low dosages seem also to be too 
high to fit the growth curve or probability integral curve. However, if the true form of 
the curve fitting such data were known, the values of the 50% end-point and of its 
standard error might be somewhat different from those found by fitting either of these 
curves to the last six points. 

3 It might be better to calculate the values of P, from the analytic expression even for 
the first arithmetic approximation to improve the value graphically determined, and 
calculation would certainly be necessary for a second arithmetic approximation. 
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4 If the arithmetic approximations were continued we should find y = 1.7787, a = 
20.53. The value of x? would be such as to give P = 0.92 for the fit to the six points, 
which would indicate that the variance of the observed from the fitted points was only 
0.29 of the amount which on the average would be due to chance or that the variation of 
the observed from the fitted points would be only 0.54 of the chance amount. 

5 For the actual computation of such a table one would not go back to the general and 
approximative equations which one would use to determine a and y for given values of 
Aand B. As the values of A and B are the same for the observated data and for the 


ce 


values fitted to them, one would start with assumed values of ‘‘y’” and ‘“‘a’’ and deter- 
mine A and B from the equations. 


A = 2m; tanh “a” (i — “‘y’”), 2B = Yin; tanh ‘‘a” (i — “y’”’). 


‘oe 


This would give a table for A and B in terms of ‘“‘y’”’ and ‘‘a’’ from which the table for 
“‘y” and ‘‘a’’ could be obtained by inversion. 











